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Abstract 

Heavy-light mesons are described in an effective quark theory with a two-body 
vector-type interaction. The bilocal interaction is taken to be instantaneous in 
the rest frame of the bound state, but formulated covariantly through the use 
of a boost vector. The chiral symmetry of the light flavor is broken sponta- 
neously at mean field level. The framework for our discussion of bound states is 
the effective bilocal meson action obtained by bosonization of the quark theory. 
Mesons are described by 3-dimensional wave functions satisfying Salpeter equa- 
tions, which exhibit both Goldstone solutions in the chiral limit and heavy-quark 
symmetry for mg — > 00. We present numerical solutions for pseudoscalar D- 
and i?-mesons. Heavy-light meson spectra and decay constants are seen to be 
sensitive to the description of chiral symmetry breaking (dynamically generated 
vs. constant quark mass). 



"Permanent address: Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear 
Research, Dubna, Head Post Office, P.O. Box 79, 101000 Moscow, Russia 
^e-mail address: yura@darss.mpg.uni-rostock.de 

'^Supported by a research fellowship of the Deutsche Forschungsgemeinschaft (DFG) 
'^e-mail address: weiss@mailserv.zdv.uni-tuebingen.de 



1 Introduction 



One of the most interesting problems of hadronic physics is the study of hadrons 
consisting of heavy and hght quarks. For example, 5-meson decays play an important 
role in determining the elements of the Kobayashi-Maskawa matrix, including the CP- 
violating phase Furthermore, rare decays of heavy mesons may indicate deviations 
from the standard model. For the description of physical (hadronic) decay processes 
one needs to know the meson wave functions and form factors. 

Heavy-light systems are also challenging from the point of view of strong interac- 
tions. The physics of light flavors is largely governed by the chiral symmetry of QCD 
and its spontaneous breaking in the vacuum. In particular, this entails the existence 
of light pseudoscalar Goldstone bosons, the pions. A simple model exemplifying the 
dynamical mechanism of chiral symmetry breaking is the Nambu-Jona-Lasinio model, 
which is based on a local two-body interaction between quarks [0, ^ ||. For heavy 
flavors (c, b) the situation is quite different. In this realm chiral symmetry breaking 
plays a minor role. The strong dynamics of heavy quarks simplifies because of the 
fact that they behave essentially like classical particles, i.e., their off-shellness in a 
bound state is small compared to the quark mass. Thus, heavy quarkonia (cc, bb) are 
successfully described by non-relativistic potential models using a Coulomb-type in- 
teraction at short distances and a linear confinement potential at large distances P, ^ . 
Heavy-light hadrons occupy an intermediate position. In heavy-light bound states, 
the velocity of the heavy quark is conserved in the limit rriQ oo, and the mass 
spectrum becomes independent of the heavy-quark spin. Recently, a general effective 
theory for heavy-light mesons has been constructed on the basis of both the chiral 
symmetry of the light flavors and the heavy-quark limit for the heavy flavors 0, ^. 
However, the coefficients of this effective Lagrangian are not determined by the sym- 
metries, but by the details of the underlying dynamics, QCD. Thus, for a quantitative 
understanding of heavy-light systems it is necessary to consider a suitable approx- 
imation to QCD at quark level. This approximation must take into account certain 
qualitative features of QCD, most importantly quark confinement and the spontaneous 
breaking of chiral symmetry. Moreover, it should incorporate the Coulomb interaction 
between quarks at short distances. These requirements make it necessary to consider 
models with non-local effective interactions, which can simulate both short-distance 
and non-perturbative long-distance effects. 

Here we consider the possibility of describing heavy-light mesons in the framework 
of a bilocal effective quark theory [T^. To incorporate the spontaneous breaking of 
chiral symmetry, we use a Lorentz- vector effective interaction motivated by QCD. A 
special feature of our approach is that we take the interaction to be instantaneous in 
the rest frame of the meson bound state []1T], |T2|, There are a number of reasons for 
this choice. Such an interaction leads to a transparent description of mesons in terms of 
Salpeter wave functions, which satisfy simple Schrodinger-type equations. Neverthe- 
less, the relativistic kinematics as well as effects of spontaneous chiral symmetry break- 
ing are incorporated exactly in the form of Foldy-Wouthuysen factors arising in the 
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reduction of the Bethe-Salpeter equation. Such a model reduces in the non-relativistic 
hmit to the potential model for heavy quarkonia (QQ)- Furthermore, instantaneous 
interactions in the form of the Coulomb gauge have long been used to study chiral 
symmetry breaking in a quasiparticle picture of the QCD vacuum [11^ [TO, [T] 



In general, potential models do not possess full relativistic invariance because the 
Fock space is restricted to gg-pairs and an instantaneous interaction is assumed. Nev- 
ertheless, in our model the interaction is written in a relativistically covariant form 
through the use of a boost vector proportional to the bound-state total momentum. 
Physically, this corresponds to the intuitive picture of a potential moving together with 
the bound state, as is used to describe moving atoms in QED. Moreover, for heavy- 
light bound states this boost vector coincides with the 4-velocity of the heavy quark, 
which provides a natural connection of this approach to heavy-quark effective theory. 
This relativistic formulation allows us to define Bethe-Salpeter amplitudes and wave 
functions for moving particles, which is crucial for the calculation of formfactors or 
decay matrix elements leading to the Isgur-Wise functions. A description of heavy 
meson and baryon weak decays in the heavy quark limit based on Bethe-Salpeter wave 
functions has been developed by Hussain et al. [ITS. Recently, Dai et al. have con- 



sidered a covariant instantaneous interaction in this context |jT9[. In contrast to these 
approaches we do not take the heavy-quark limit from the start but formulate the 
description of bound states for arbitrary quark masses, discussing the simplifications 
arising in the heavy quark limit afterwards. An advantage of this covariant approach is 
precisely the possibility to describe in a unified manner both light and heavy mesons, 
using appropriate potentials. For example, for light quarks one may employ a separable 
approximation to an intermediate-range potential in momentum space. In this case 
one obtains a consistently regularized version of the Nambu-Jona-Lasinio model, in 
which the 3-dimensional momentum space cutoff moves along with the bound state 



20|| . We should point out that the relativistic formulation using the boost vector is 



useful also in the light quark sector, e.g. in defining the pion decay constant pT| . 

Nowak et al. have established the general relation between dynamics at quark level 
and heavy quark effective theory by performing a gradient expansion of the fermion 
determinant in the heavy quark limit [^. The bilocal meson action derived in our 
approach may serve as an explicit realization of their ideas. However, we shall consider 
here the full momentum-dependent theory and do not restrict the effective action to 
the long- wavelength limit. 

In this paper, we first develop the formal framework for the description of meson 
bound states in a bilocal effective model with a covariant instantaneous interaction, 
which takes into account the spontaneous breaking of chiral symmetry. We then apply 
this model to the study of light-light and heavy-light pseudoscalar mesons. Specifically, 
we want to demonstrate that the momentum dependence of the constituent quark 
mass, which is a consequence of the dynamical nature of chiral symmetry breaking, 
has important effects on the spectrum and decay constants of heavy-light mesons. 
The description of meson transitions and form factors in this unified approach will be 
given in the future . 
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This paper is organized as follows. In sect. 1 we formulate the effective quark the- 
ory and introduce the instantaneous interaction in the covariant formulation |jTT|, |T^ . 
We discuss the phenomenological 3-dimensional potentials considered in the appli- 
cations. We then bosonize the model and obtain an effective bilocal meson theory. 
The Schwinger-Dyson equation for the vacuum and the Bethe-Salpeter equation for 
the meson fluctuations are derived from the effective meson action. We outline the 
calculation of matrix elements between bound states needed to describe meson decay 
constants or semileptonic decays. In these derivations we assume a general covariant 
interaction kernel. In sect. 2 we then discuss in detail the equations describing bound 
states for the case of a covariant instantaneous interaction. The Schwinger-Dyson 
equation describes the quark quasiparticle spectrum. The Bethe-Salpeter equation for 
mesons is rewritten in terms of a Salpeter wave function. We make use of a Foldy- 
Wouthuysen transformation to simplify the form of the resulting equations. In sect. 
3, we then apply the model to the description of light-light and heavy-light mesons. 
We discuss the chiral limit, — >• 0, and study the behavior of the pseudoscalar me- 
son in this regime. We then investigate heavy-light bound states {Qq}- In the heavy 
quark limit, mq oo, reduced bound state equations are obtained, which exhibit the 
heavy-quark spin symmetry. We present numerical solutions of the full equations for 
pseudoscalar D and _B-mesons. Specifically, we compare meson properties calculated 
with the momentum-dependent quark mass from the Schwinger-Dy son-equation with 
those obtained using a constant light quark mass. We find that heavy-light meson 
properties are sensitive to whether chiral symmetry breaking for the light flavor is de- 
scribed self-consistently or by a constant constituent quark mass. A summary and an 
outlook are given in sect. 4. 

Appendix A.l deals with the partial-wave analysis of the meson Salpeter equation. 
In appendix A. 2 we describe the numerical solution of the partial-wave equations using 
the Multhopp method, which has been very sucessful in the context of constituent quark 
models |^^. In particular, we consider it necessary to comment on the complications 
presented by the momentum-dependent Foldy-Wouthuysen factors in the bound state 
equations of our model. Appendix B outlines the evaluation of the matrix elements 
for the normalization of meson wave functions and the calculation of meson decay 
constants from the bilocal effective meson action. 

2 Effective quark theory with instantaneous inter- 
action 

The basis of our description is a quark theory with an effective two-body interaction 
in the color-octet channel. It is defined by the action of the quark field, 

W = Jd'xq{x)G,\x)q{x) - JJd'xd'yj'^^{x)Df,{x-y)f\y). (1) 
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Here, Go is the free quark Green function, 

= i^-m°, (2) 

where mP = diag(m5', . . . , mPp^^.) is the quark current mass matrix, with Nf the number 
of flavors. The quark color current is given as 

fJx) = q{x)(^)^,qix) , (3) 



2 

where A*^ are the Gell-Mann matrices of SU{3)c. The bilocal interaction kernel, 

K(^-y) = 5'^V^(^-?/), (4) 

can be thought of as an effective gluon propagator, which describes part of the non- 
abelian effects of QCD in a phenomenological way. Note that the interaction is of 
Lorentz-vector type, as motivated by QCD, and thus chirally invariant. 

For the purpose of describing meson bound states we rewrite the interaction term 
of eq.(0) in the form 



d^xd'^y qB{y)qA{x)kAB,CD{x - y)qD{x)qc{y), (5) 



with the kernel 
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KabM^-v) = li^uil,)tsj:^^^iihj^-^D{x-y). (6) 

a=l 

Here, A = {r,a,i}, . . . , D = {u, 6, 1} are a short-hand notation for the Dirac spinor, 
color and and flavor indices. In the following we want to study meson [qq) bound 
states. We therefore make the well-known color Fierz rearrangement |p5| 



8 4 2 3 

a=l p=l 

where eap^y is the antisymmetric Levi-Civita tensor. This identity allows to rewrite 
the interaction completely into the attractive color-singlet [qq] and antitriplet {qq) 
channels, while the repulsive color-octet and sextet channels are absent in a natural 
way. We consider only the color-singlet part of eq.(|^). The relevant part of the action, 
eq.(|l|), can then be represented in the form 

W = JJd'xdS!^{q{y)q{x)){-G^\xMx-y) (8) 
1 



{q{y)q{x)) K{x - y) {q{x)q{y)) 

with the color-singlet interaction kernel 

K{x-y) = -f^(^-f^g^D{x-y). (9) 
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Here, the bilinear {q{y)q{x)) is contracted over color indices, and A^^c = 3 is the number 
of colors. 

If the model action eq.(|l]) is considered as a Euclidean field theory with a covariant 
gluon propagator, eq.(|l]) is known as the so-called Global Color Model In that 



approach, bound states are studied through the Euclidean correlation functions. Sim- 
ilar covariant models have been investigated in [^, Here, we consider eq.(|I|) in a 
different context. For reasons already stated above we wish to have an effective theory 
with an instantaneous interaction. Such an interaction leads to a simple description of 
bound states in terms of 3-dimensional wave functions in Minkowski space, which for 
heavy-heavy bound states naturally reduces to the successful non-relativistic potential 
model. The concept of an instantaneous interaction can be formulated covariantly by 
letting the potential move along with the bound state. 

Consider a bound state of two quarks interacting through a bilocal effective in- 
teraction of the form eq.(P). From the principle of translational invariance it follows 
that one can separate the center-of-mass motion of the bound state from the relative 
motion in the form of a plane wave. The momentum of the c.o.m. motion equals the 
total momentum of the quark pair, V^. Given this, one can more or less arbitrarily 
define a conjugate coordinate, X^, representing the absolute position in space-time 
and a relative coordinate, z^, to describe the internal structure of the bound state. 
The condition is that under a translation by a constant vector a, x — x + a, ?/ — ?/ + a, 
these coordinates transform as X — ^ X + a, z — >■ ^. This is satisfied for any linear 
combination X = ax + {1 — a)y and z = x — y. We shall take a = | in the following^. 

We now want to substitute the general bilocal interaction kernel of eq.(||) by an 
instantaneous one, i.e., by a potential. This can be done in a relativistically covariant 
way by replacing eq.(|^) by 



K'^ix^y) = K'^l^x-y 

Here, 



fl^flV{z^)6{z\\). (10) 



4 = Vf^i^-V), = Zf,-z^^, (11) 

and r]^ is a boost vector proportional to the total momentum eigenvector of the bound 
state, Vfj,, 

V 



In eq.(|I^), the 5-function, 6{z-ri), guarantees the instantaneousness of the exchange 
interaction in the rest frame of the bound state. The transversality of the exchange 
interaction is ensured by the fact that the 3-dimensional potential is a function only 

""^The definition of the center-of-mass coordinate a.s X = [raix -\- m2y) / {mi +m2) has significance 
only in the non-relativistic hmit. In a relativistic theory it is in general impossible to write the 
two-body hamiltonian as the sum of two terms describing c.o.m. and relative motion. 
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of the perpendicular part of the relative coordinate, V{z-^). The sign in eq.(|r^) has 
been chosen such that an attractive 3-dimensional potential leads to to an attractive 
interaction in the color-singlet (and antitriplet) channel. Eg . (p!0D describes a potential 
moving together with the bound state. This form of interaction leads to bound-state 
amplitudes with well-defined Lorentz transformation properties. We remark that one 
can arrive at the form eq.(|TU]) by studying moving bound states within a canonical 
quantization approach to gauge theories 

For the bound state at rest one has V = 0, so that r^^ = (1, 0). In this frame the 
kernel takes the form 

ir(i'°)(a;,y) = -7o®7o^(x-y)5(xo-2/o). (13) 

In electrodynamics this kernel corresponds to the usual Coulomb gauge, with V{:x. — y) 
the Coulomb potential. The form eq.(|13D has been widely used as a model to describe 



the breaking of chiral symmetry in strong interactions |T^, [16|, |T7 . 

Let us now assume that the bound state contains a heavy quark (c, b) and a light 
antiquark [u, d, s). In the limit of heavy quark effective theory, mq —>■ oo, one has 

= ^ ^ (14) 

where v is the 4- velocity of the heavy quark, f ^ = 1. Thus, in this limit the interaction 
kernel, eq. (p!OD , takes the form 

K'^ix-y) = -i}(^ijV{z^)6{v-z), (15) 

where = z — v{v-z). The fact that the interaction kernel explicitly involves the heavy 
quark 4-velocity leads to a natural relation of this potential approach with heavy-quark 
effective theory 0, |^ . 

The transverse potential, V{z-^), is a phenomenological input to this model. Its 
form may be chosen depending on the type of bound state one wishes to study (light- 
light, heavy-light). We shall in the following work with a potential, which in the rest 
frame is of the form 

4 (T 

\/(x-y) = --^ + a\ (16) 
6 r 

with r = |x — y|. Note that when expressing the interaction in an arbitrary frame one 
must take into account that z-^ is subject to a non-Euclidean metric. The first term of 
eq.(|T^ is the Coulomb potential describing one-gluon exchange. In the investigations 
in this paper we shall take as as constant. One could also employ in this approach an 
asymptotically free potential with as = a^dq^) |jl6|, [l^. The linear potential in eq.(|l6D 



implements quark confinement in a phenomenological way. In our study of heavy-light 
mesons we use parameters determined within the non-relativistic potential model for 
charmonium, as ~ 0.3 . . . 0.5, a ~ 0.4 GeV ^. For reference, we shall also consider 
the oscillator model of Le Yaouanc et ai, V(x — y) = VqV^, with Vq ~ 0.2 . . . 0.4 GeV, 
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which gives a good overall fit to the light meson mass spectrum |]2T|. We remark that 
one may include in eq.([T6|) also an intermediate-range attractive potential, which can 
be thought of as a crude representation of instanton effects |3^, ^ . Another possibility 
is a separable potential in momentum space, which leads to a Nambu-Jona-Lasinio- 
type model. Such an interaction has been considered within this covariant approach 
in m. 



3 Bound states 



3.1 Hadronization 



For the discussion of meson bound states it is convenient to "hadronize" the quark 
theory, i.e., to formally rewrite it as an effective meson theory p, |1^. For the general 
discussion here we assume for simplicity a general covariant interaction kernel. The 
equations obtained will then be specified to the case of a co-moving instantaneous 



interaction, eg . ([T0|) , in the following two sections. 



Let us consider the functional integral 

Z = J VqVq expiW[q,q]. (17) 

Here, W is the color-singlet part of the effective quark action, eq.(D, which in symbolic 
notation can be written as 

W[q,q] = (qq,-G^^) + ^{qq,Kqq). (18) 

After integrating over the quark fields with the help of the Legendre transform one 
obtains 

Z = JvMexptWeslM], (19) 
with the effective meson action 

W,s[M] = N,l^-l{M,K-'M)-iTT\og{-G,'+M)y (20) 

Here, A4 = A4ij{x,y) ~ (li{x)qj{y) is a bilocal meson (color-singlet) field. It has the 
structure of a matrix in Dirac spinor and flavor space. In eq. (pO[) , the symbol Tr implies 
integration over the continuous variables as well as the traces over spinor and flavor 
indices. 

The vacuum of the effective meson theory is given as the minimum of the effective 
action, eq.(^). The condition of minimum reads 
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The vacuum solution of this equation is translationally invariant. Let us denote it by 
(S — m*^). Then we obtain from eq.(pID the Schwinger-Dyson equation 

S = m'^ + iKG^, (22) 

where 

G^^ = ip-E. (23) 

Mesons are described as fluctuations of Ai around the vacuum configuration. Ex- 
panding the action, eg . (|20|) , around the minimum, with = (S — m") + Ai', one 
obtains 

W,s[M] = Wesin 

+ nJ --{M', K-^M') - -Tr (Gs-M')' - ^ E -Tr {-Gj^M'T ■ (24) 
L 2 z n ) 

The vanishing of the second variation of this effective action with respect to 



eff 



■ r = 

M'=0 



6M'6M' 

leads to the homogeneous Bethe-Salpeter equation for the vertex function, F, 

r = -tK{Gj:TGj:). (25) 

This equation describes the bound state spectrum. It corresponds to the usual Bethe- 
Salpeter equation in ladder approximation. 

Given the solutions of eqs.(^, ^) describing the bound state spectrum, one can 
calculate matrix elements between on-shell bound states describing meson decays, semi- 
leptonic processes, etc., in the framework of the effective meson action, eq. (pOD . The 
bosonized action summarizes these processes in a concise way and greatly facilitates 
the evaluation of the corresponding matrix elements. One formally expands the bilocal 
field, J\4, in bound state amplitudes, 

x + y\ ^ f dV^ 



M{x,y) = M(x-y ^) =J2j 



(27r)3/2V2^ 

X / ^e-(-^){e^^^a^,(P)r^(g|P) + e-^^'-^ au{V)T n{q\V)]. (26) 

The sum runs over the set of quantum numbers, if, of hadrons contributing to the 
bilocal field. Here, the individual bound states have mass Mh and total 4-momentum 
V = (UH^V), with luhCP) = (V^ + Mlf'^. The amplitudes r^^(g|P) and TH{q\V) = 
^H{q\ — V) are on-shell solutions of the Bethe-Salpeter equation, eq. (|25|) , with 3- 
momentum V. The coefficients a^{V),aH{V) may then be interpreted as creation 
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and annihilation operators of physical (on-shell) mesons, and matrix elements can be 
calculated as usual. In particular, the bound states amplitudes are normalized by the 
requirement that the matrix element of the quadratic (free) part of the effective action, 
eq.®, 

lyf = -i\Nj:i{G,,Mf (27) 

have the normalization corresponding to a physical (elementary) particle. This is sim- 
ply the statement of the correct relativistic dispersion law of the c.o.m. motion of the 
bound state0. 

In order to describe meson weak decays, we couple the theory eq. (|l|) to an external 
weak leptonic current. The corresponding lagrangian at quark level is given by 

>Csemi = ^{V,,{Qi{x)0^qj{x))l^{x) + h.c.}. (28) 

Here, is the leptonic current, 

l^{x) = l{x)Of,ui{x), I = (e,/i,r), = (z/e, z/^, z/^), (29) 

= 7m(1 + 75), 

Vij are the elements of the Kobayashi-Maskawa matrix, and Gp is the Fermi constant. 
In eq. (|28|) , Q denotes the column of (m, c, t)-quarks, q the {d, s, 6)-quarks. On the level 
of the effective meson action, eq. (|20|) , the electroweak coupling can be incorporated by 
shifting the bilocal field by the local leptonic current, 

Mij{x,y) Mij{x,y) + Lij{x,y), (30) 
L,^{x,y) = 9Zs\x-y)Kp\{x)e'''-'^, (31) 

where Vl is the momentum of the leptonic pair. After this shift, the part of eq.(P7[) 
describing meson decay into leptons is 

W2L = -^iV,Tr(Gs-MGsL). (32) 

Semileptonic decays of mesons, which for heavy-light mesons are parametrized in terms 
of the Isgur-Wise functions are mediated by 

W^iemi = tN,Ti{G^MG^MG^L). (33) 



In the following we consider in detail the Schwinger-Dyson-equation, eqs. (^2D , 
and the Bethe-Salpeter-equation, eq. (P^D , for a covariant instantaneous interaction, 
eq.(|10D. In particular, bound states will be described by 3-dimensional wave func- 
tions. The evaluation of meson observables then proceeds in general as follows. Given 

^For a covariant instantaneous interaction, eq.(|l^), the interaction does not contribute to the 
normahzation of the bound state, cf. appendix B. 
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the solution for the bound state wave function in the rest frame, one can reconstruct 



the Bethe-Salpeter amphtude in the rest frame, cf. eqs. (|^, |58|, jG^) below. By virtue 
of the relativistic formulation of the interaction, eg . (p!0D , the bound state amplitude 
transforms covariantly and may be boosted to an arbitrary velocity. Matrix elements 
of the vertices eqs.(|3^, |33D etc. between moving bound states can then be calculated 
from the expansion of the bilocal field, eg. . The evaluation of the matrix element 



for pseudoscalar meson decay in this approach is discussed in appendix B. 
3.2 The quark spectrum within a meson 

For an instantaneous interaction of the form eq.(p!0|) the Schwinger-Dyson equation. 



eq.(P4), describes the spectrum of quark quasiparticles in the rest frame of the meson. 



In momentum space it takes the form 

S(p^) = rh'-^J^,Vip^-q^)f,GMfl■ (34) 



Here, the 4-dimensional Fourier transform of the interaction kernel, eq.(|TO|), depends 
only on the transverse part of the relative momenta. 



V{p^-q^) = J d''xe-'^P-''^''V{x^)5{x-r]). (35) 
In the following we consider eq.(|3^ in the rest frame, where = (1,0), fj = 7°, and 



eq.(|35D coincides with the usual 3-dimensional Fourier transform. Assuming that S(p) 
is diagonal in flavor, S = diag (Si, . . . S^r^), eq.(|5^) splits into identical equations for 
the S„ with bare mass m°, n = 1, . . . Nf. We shall omit the flavor index on S in the 
following. The quark self-energy in the rest frame has a scalar and a vector part, 

S(p) = A(|p|)|p| + i?(|p|)p-7. (36) 

In the following it will be convenient to make a polar decomposition of the quark energy 
in the form 

S(p) + p-7 = E{\p\)S\p), (37) 

where 5'^(p) is the square of a Foldy-Wouthuysen matrix, 

S^\p) = sin0(|p|) ± p-7cos0(|p|), (38) 
E(|p|)sin0(|p|) = A(|p|)|p|, E(|p|)cos0(|p|) = (l + 5(|p|))|p|. (39) 

Here, < 0(|p|) < |vr is the so-called chiral angle. (In the quasiparticle language this 
angle deflnes the rotation of the massive quasiparticle spinors relative to the free quark 
spinors in the vacuum of broken chiral symmetry [|l^, |15| . ) In this parametrization the 
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quark propagator, Gs{q), becomes 

GM = = ( t'l ^ + A-'t . h» (40) 

^-S(q) \qo-E{q)+%e qo + E{q) - le/ 



7o 



A+(q) , A_(q) 



.go-^(q) + '^e qo + E{q)-ie 
Here, the matrices A±(q), A±(q) are defined as 



where 



A±(q) = S-\q) L S{q), A±(q) = S{q) A± ^-^(q), (41) 



A± = i(l±7o) (42) 



is the usual projector on positive and negative energy components. The first power of 
the Foldy-Wouthuysen matrices, S'^^(q), can be expressed from eg. (|38D 



as 



S^\q) = cosz/(|q|) ± q-7sinz/(|q|), i^(|q|) = |(|7r - 0(|q|)). (43) 

Inserting the quark propagator in the form of eq. (^OD into the Schwinger-Dyson equa- 
tion, eq. (|34D , and taking traces one obtains a system of equations for -E(|p|) and 0(|p|), 



(44) 



^^^y(p-q) sin0(|q 



E(|p|)cos0(|p|) = IpI - l|^y(p-q)p.qcos0(|q|). 



{2ny 

These equations define the single-particle spectrum of two quarks forming a bound 
state. From eqs.(H3) one obtains an integral equation for the chiral angle. 



m° cos 0( IpI) — IpI sin 0( IpI) 
1 f d^q 



2 



i2ny 



V{p-q) [cos0(|p|)sin0(|q|) - p-q sin0(|p|)cos0(|q|)] . (45) 



Here, 0(|p|) satisfies the boundary conditions 0(0) = ^vr and 0(|p|) — > for |p| oo. 
From the chiral angle, the quark energy can be obtained as 

E{\p\) = m°sin0(|p|) + |p|cos0(|p|) 

1 f d^q 



2 



i2ny 



V{p-q) [sin0(|p|)sin0(|q|) + p-qcos0(|p|)cos0(|q|)]. (46) 



We remark that absolute values of the quark energy may be shifted by a constant 
amount, Eq, without affecting the bound state spectrum, by adding to the potential a 
contact term, Vo(p — q) = — -E'o(2vr)^5*^^)(p — q), cf. the discussion in [1^, |T^. 
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If the potential contains a Coulomb interaction, Vc(p — q) = — |as/(p — q)^, the in- 
tegrals in eqs. (^4]) are ultraviolet divergent and require renormalization. Consequently, 
we perform in eqs.(^4D a wave function and mass renormalization and replace the bare 



quark kinetic energy and mass, |p| and m , by 



iPl 



m 



iPl 



m 



(27r)3 



m 



qP + m' 



02 



(47) 
(48) 



The equations resulting from eqs. (^ , ^q) upon this substitution are finite. The 
boundary conditions on the chiral angle remain unchanged. Furthermore, with the 
definitions eqs.(|4^, ^ ) the single-particle energy for |p| ^ oo behaves like that of a free 
quark. In the case of an asymptotically free potential, = a^dpl) it has been shown 
that the Schwinger-Dyson equation can be renormalized consistently by requiring that 
the quark propagator reduce to the free propagator at some large 3-momentum Pren, 
with |pren| — * oo . Here, we are working with a potential, the UV-divergent 

part of which is of pure Coulomb form, as = const., so that the definition eq.(^7D is 
sufficient^. The linear potential in eq.(|T^) is UV-finite and does not contribute to the 
renormalization. The redefinition of m°, eq . (^8]) , is not a renormalization in the strict 
sense, as it corresponds to using a momentum-dependent renormalization constant. In 
the framework of this potential model we take the point of view that eq. (|48|) is simply 



a prescription to softly break chiral symmetry for the light flavors with parameter 
m°. Hence the value of is not inferred from the QCD current mass but may be 
determined phenomenologically e.g. by fitting the pion mass. This is usual practice 
in the framework of the Nambu-Jona-Lasinio model 0, |^. Definition eq.(|^) leads to 
a smooth behavior of the pion mass and decay constant in the chiral limit, as can be 
seen from figs. 4 and 5. For heavy quarks, we will use the approximation of a constant 
constituent quark mass, cf. below. 

The ratio of the scalar to vector part of the quark self-energy. 



m(|p|) 



^(IpI)IpI 

l + Bi\p\) 



|p|tan0(|p|). 



(49) 



can be interpreted as a momentum-dependent "constituent" quark mass, which reduces 
to the current quark mass, m°, 

Eq.(H^) constitutes a non-linear integral equation for the chiral angle. 



in the limit |p| ^ oo ||17|| . 

IpD- By 

performing the angular integral one obtains a 1-dimensional equation, which we solve 
numerically using the relaxation methods of Adler et al. |TB| 

Numerical solutions of the Schwinger-Dyson equation for the chiral angle, 0(|p|), 
and the quark energy, i?(|p|), for a light quark fiavor are shown in figs. 1 and 2. Here, 

^For a pure Coulomb potential without a mass scale the definition eq.(^) falls within the scope 
of multiplicative renormalization. In other words, for a pure Coulomb potential the subtraction 
procedure of Finger and Mandula flj] coincides with the renormalization of Adler and Davis . 
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the potential is the one used below to describe D and 5-mesons (a = 0.41 GeV, as = 
0.39). Note that the attractive interaction causes the quark energy to become negative 
at small momenta. This fact is intimately related to the spontaneous breaking of chiral 
symmetry and the emergence of Goldstone bosons [|T5[. The dynamical quark mass, 
m(|p|), is shown in fig.3. We also show the corresponding quantities for the oscillator 
model of []T3| (Vq = 0.247 GeV), for which eq.(^3|) reduces to a differential equation for 

0(|P|)- 

For heavy quark flavors, the bare quark mass, m°, is much larger than the dynam- 
ically generated mass. One may then neglect the dynamical contribution to the mass 
and approximate S(q) by a constant constituent quark mass, S(p) = mq. In this case 



we have in eq.(|39D A(|p|)|p| = mg, -B(|p|) = 0, so that 

0(|p|) = arctan^, E{\p\) = J\p\^ + ml, u{\p\) = ^arctanM. (50) 



The chiral angle and the quark energy for a heavy quark are also shown in figs. 1 and 
2. A constant quark mass is also frequently used to describe the light-quark sector, 
see e.g. We will comment on the consequences of this approximation below. 

In the heavy-quark limit, if mg becomes much larger than the momentum of 
the relative motion of the heavy-light bound state, one may approximate the Foldy- 
Wouthuysen angle of the heavy flavor by its value at |p| = 0, 0(|p|) = ^tt, z^(|p|) = 0. 
Furthermore, the heavy-quark energy becomes E{p) ~ mq + ^|pp/mQ ~ niq + 
0{l/mq)] to leading order one can neglect the heavy-quark kinetic energy. In other 
words, in this approach the heavy-quark limit amounts to neglecting the variations of 
the heavy-quark energy and chiral angle over the range of momenta contributing to 
the bound-state wave function. 



3.3 Equations for bound state wave functions 

Given the quark single-particle spectrum in the vacuum of the effective theory, we now 
turn to the description of meson bound states. For the covariantly written potential 
kernel, eq.(|lOl), the Bethe-Salpeter amplitude, T{p\V), depends only on the transverse 
part of the internal momentum, p-^. With this interaction the Bethe-Salpeter equation, 
eq.(p5D, for a bound state of flavor qiqj in an arbitrary frame reads 

nP^\V) = ^ j ^V{p^ -q^)iG,{q+\V)T{q^\V)GM-\m- (51) 

Here, Gn = Gj^^ is the quark propagator of flavor n defined in eq. (|40|) , and V de- 
notes the total momentum of the bound state. For an instantaneous interaction it is 
convenient to describe the bound state in terms of a Salpeter wave function, 

^{q^\V) = ^j'^GM+\V)mV)GM-\n (52) 
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where q^^ = q-r]. In the following we shall restrict ourselves to the rest frame, where 
V = (M, 0), with M the bound state mass, and = go- We drop the argument V in 
r(p|P), \l/(p|'P) in the following. Performing in eq. (|5T|) the integral over go one obtains 
for ^^(p) the equation 

Mr( ^ n+_(p) n_+(p) 

Here, Ep = Ei{\p\) + Ej{\p\) is the sum of the single-particle energies for the two quark 
flavors, and 

n±^(p) = A.±(p)7or(p)7oAiT(p), (54) 
r(p) = /^1^(P - q)7o*(q)7o, (55) 

with the projectors An±{p), A„^(p) for flavor n defined in eq.(5l). Rather than working 
with eq.(p3D directly it is more suitable to introduce an "undressed" wave function, 

&(p) = 5,(p)vl/(p)5,(p). (56) 

Here, Sn{p) are the Foldy-Wouthuysen matrices, eq.(^3|), corresponding to the two 
quark flavors. The new wave function satisfies an equation analogous to eq. (^3]) , 



I^(P) = ni_M + n^ 

E„-M E„ + M' ^ ^ 





where H±^ (p) now involves the free projectors, eq.(|4 

ri±^(p) = -Lsr\p)np)Sr\p)A^. (5t 
As a consequence, the new wave function satisfies the constraints 





A±^ (p) A± = 0. (59) 



By applying projection operators, eq. (|57D can be written in the form of a Schrodinger 
equation, 

[Ep T M] A±|r (p) = n±^ (p) . (60) 

It has the structure of a generalized eigenvalue equation for the bound state wave 
function, v]> (p), with M as eigenvalueQ We now decompose ^ (p) as 

(1) (2) 

^ (P) = * (P) + 70- ^ (P)- (61) 



"^The mathematical structure of the instantaneous Bethe-Salpeter equation has recently been an- 
alyzed in |3i,|3|]. 
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Because of the constraint, eg. , we can expand the functions \1> (p) completely in 
the set of Dirac matrices {75, ^7}. We write 

c^) 

^ (p) = L«(p)75 + zN('=)(p)-7 (fc = l,2). (62) 

By taking traces of eg. (|60|) one obtains after a straightforward but lengthy calculation 
a system of eguations for the component functions, L'^'^^(p), N'^'^)(p), k = 1,2. In 
calculating the traces it is convenient to introduce a dreibein in momentum space, 
{p, Cp, Cp}, with p-Cp = 0, Cp-Cp = 6°''' {a,b = 1,2) and expand the wave function in 
the set {75, iCp-'j, ip-^}- We find 



{2\ n\ 



ML''\p) - E,L''\p) = /(^i^P-q) (63) 



X 



{[cM + ^J4p-q]^ (q) + 4<p-(qxN ^ (q))} 



{2\ n 



MN ' (p) - E,N ^ (p) = / 7^1^(P - q) (64) 



X 



(2^ 

\ L^p p q ~ -'p-'q ~ ^p^q ~ ^p ^q -'p^q 

(1) 

-sjsjpx (qx .) + sJs±p(q- . )] N ' (q) 



pxqL^ (q)} 

Here, , P^ are the transverse and longitudinal 3-dimensional projectors, 

P^ = p®p, Pp^ = 1-Pp^ = -px (px .), (65) 

with a similar definition for q. The factors , are the result of the Foldy-Wout- 
huysen transformation, eg. (pB|) . We have introduced the short-hand notation 

Sp = sin z/jp cos Ujp ± cos z/jp sin Vjp = sin(z/jp ± Ujpj, 

with a similar definition for c^,s^. The angle z/„p = z/„(|p|) for flavor n = i,j has 
been defined in eg. PB|) . The functions Ep = Ei{\p\) + Ej{\p\) and Uip, Ujp contain the 
entire information on the single-particle spectrum. They have to be provided either 
as solutions of the Schwinger-Dyson eguation for the guark self energy, eg. (|4^) , or by 
the approximation of a constant guark mass, eg. (pOD . 

The eguations egs. (|63| , |64D for the bound state wave functions have a very compact 
form. Note in particular that the pseudoscalar-axial and scalar-vector meson mixing 
is taken into account here. Corresponding eguations have been given by Le Yaouanc 



Cp = cos Vip cos Vjp =F sin z/jp sin vjp = cos(z/jp ± Ujp), (66) 
± 
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et al. for the case of equal quark masses and a harmonic oscillator potential |2T| . 



Through the use of the Foldy-Wouthuysen transformation, eq. (|56|) , we have preserved 
the simple structure of the wave function, eg . (|59D , even in the general case of unequal 
quark flavors. Note that in the case of identical quark flavors (isospin limit), one 
has s~ = in eq.(|66D, and the equations for the L- and N-component decouple]^. 
We remark that equations in a parametrization somewhat different from to that of 
eqs.(|63|, |64D have been derived by Lagae using a variational approach The use of 
the Foldy-Wouthuysen transformation in the form eq.(|5BD provides a simple and more 
transparent alternative. 

The bound state wave functions are normalized by condition eq.( [B.8| ), which fixes 
the relativistic dispersion law for the bound state. From eqs. (p6|, p^ , |60|) the Bethe- 
Salpeter amplitude can be reconstructed from the "undressed" wave function in the 
rest frame. The wave functions contain the full information about the bound state 
on its mass shell. In appendix B, the pseudoscalar meson decay matrix element is 
evaluated in terms of the wave function. 

For calculation of the meson spectrum it is necessary to decompose eqs.(|63|, |6^) 
into equations for bound states of given angular momentum and parity in the rest 
frame. This is done in appendix A.l. The resulting radial equations are solved numer- 
ically using the Multhopp method [|24|. In the presence of the momentum-dependent 



Foldy-Wouthuysen factors in the integral equation this momentum-space method is 
much more convenient than the commonly used matrix methods, in which the matrix 



elements of the potential are evaluated in position space |34, 35|. In appendix A. 2 we 



show how the Multhopp method can be adapted so that essentially no complications 
arise from the Foldy-Wouthuysen factors. 



4 Light— light and heavy— hght mesons 

The bound state equations derived from the bilocal effective action, eq. (|20D , describe 
mesons for arbitrary quark masses. We now apply this description to the study of 
light-light and heavy-light mesons. Before embarking on the numerical solution of the 
Salpeter equations it is worthwhile to investigate some limiting cases. In particular, we 
wish to demonstrate how eqs.(|5B|, |5^ describe both pseudoscalar Goldstone bosons in 
the chiral limit, m° — >■ 0, and degenerate heavy-light mesons in the heavy-quark limit. 
Let us first consider eq. (|63D in the isospin limit with two light quark flavors, i.e., for 
the pion. In this case one has Uip = ujp = i^{\p\) = |(|7r — 0(|p|))) so that 

c+ = sin0(|p|) s+ = cos(/)(|p|) 

Cp = 1 Sp = 

with 0(|p|) the solution of eq. (^) . If furthermore = 0, comparison of the Salpeter 
equation, eq.(^), with the Schwinger-Dyson equation, eq.f^if), shows that eg . (|63D 



^The equations for wave functions written down in ||l3|| apply only to this case of equal quark 
flavors, since the L-N mixing term has been omitted there. 
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possesses a massless Goldstone mode solution of the form 

L(^)(p) oc sin0(|p|), L^^Hp) = 0, M = 0. (68) 

For this property it is crucial that chiral symmetry is broken spontaneously by the 
same interactions which bind the quarks in the meson, i.e., that the quark self-energy 
is taken as the solution of the Schwinger-Dyson equation and not approximated by a 



constant constituent quark mass ||3T|]. Note also that the pseudoscalar wave function 
is normalizable only for m° > 0. 

Of particular interest is the case of a small but finite current quark mass, the 
chiral limit. Specifically, we want to see how the pseudoscalar meson mass and decay 
constant behaves in the limit m° 0. We thus solve the Schwinger-Dyson equation, 
eq.(Q), and the pseudoscalar bound state equation, eqs. (|63|, |A.7| ), for small current 



quark masses, = m° = m°. Results for the linear plus Coulomb potential as well 
as the oscillator model are shown in figs. 4 and 5. As can be seen from fig. 4, the pion 
mass vanishes in the chiral limit like oc m°, in accordance with current algebra. The 
pion decay constant, /tt, is shown in fig. 5. As expected, it approaches a finite limit if 
m° — i> 0. We remark that for a pure oscillator potential (Vq = 0.247 GeV) we reproduce 
in the chiral limit the value of Le Yaouanc et ai, f-j^ = 20 MeV, while for a pure linear 
potential (cr = 0.40 GeV) we find /^r = 0.11 MeV, which is in agreement with the value 



obtained by Adler and Davis |16]. Pion properties in the chiral limit have also been 
studied by Alkofer and Amundsen using the inhomogeneous Bethe-Salpeter equation 
T7[] . For an instantaneous interaction, the wave function description of bound states. 



eq.(|5^, is a convenient alternative, especially for the study of excited states. Also 
shown in fig. 5 is the decay constant for the radially excited pion state, which vanishes 
in the chiral limit, as expected on general grounds The reason for this is that as 
the pion wave function approaches sin(/)(|p|), the excited state wave functions become 
orthogonal to sin(^(|p|) and thus the integral for the pion decay constant, eq. ( |B.15D , 
vanishes. The pion radial wave function, cf. eq. (|A.3|) , is shown in fig.6, for a current 
mass of m° = 1 MeV. 

It is well-known that a charmonium (linear plus Coulomb) potential with usual pa- 
rameters underestimates the strength of spontaneous chiral symmetry breaking, which 



manifests itself in the too small value of /^r [|T6l. Improved values can be obtained 



by including either an intermediate-range attractive potential PO] or transverse gluon 
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Thus, we have verified that the bilocal effective meson action, eq.(^), reproduces 
the successful phenomenology of the pion as a Goldstone boson, if the breaking of chiral 
symmetry is described consistently with the interactions which form the bound state. 

Let us now consider heavy-light mesons. Specifically, we shall investigate a bound 
state of a heavy quark and a light antiquark (Qiqj). For the heavy quark we neglect 
spontaneous chiral symmetry breaking and approximate its energy by a constant mass, 
cf. eq. (|50D . In order to take the limit mg oo, we subtract from Ep and M in eqs.([63|, 
|n^ ) the heavy quark mass, i.e., we consider binding energies relative to the heavy quark 
mass. For the light flavor, we take take into account the dynamical breaking of chiral 



18 



symmetry and take i?j(|p|), z^j-dpl) as the solution of the Schwinger-Dyson equation, 
eq.®. 

In particular, in the heavy-quark limit, if rriQ becomes much larger than the range 
of momenta contributing to the integrals over the bound-state wave function we may 
simply take i/j(|p|) = 0. In this case we have 



1 



cos(±z/jp) 
sin(±z/jp) 



cos [h'jp) = 

±sin(z/jp) = 



±S 



(69) 



P5 



while Ei{\p\) = 
equations for the 
solutions satisfy 

L«(p) 



pP = niQ + 0{l/mQ). Consequently, in eqs. (|63|, |6^ ) the 
1)- and (2)-components of the wave function coincide, and the 



L(2)(p) = L(p), N«(p) = N(2)(p) ^ N(p). 



(70) 



Eqs.( 



thus simplify to 



(M-Ep)L(p) = J^s^P-^) 



(71) 



(M - Ep)N(p) = J ^,V{p - q; 



[[cpCq + SpSqp-q\L{q) - SpSgp-(qxN(q 

W) 

:{[CpC, + SpS5((p-q) - (p X q) X . ))] N(q) + SpS,pxqL(q)} 



(72) 



The reduced equations eqs. ([7I|, |72|) exibit the heavy-quark spin symmetry, i.e., the 
pseudoscalar bound state is degenerate with the vector, the scalar with the axial vec- 
tor This symmetry is a consequence of the fact that the dynamics becomes 
independent of the heavy quark spin in the limit uiq — >■ cxo. We can demonstrate 
this explicitly for S-wave bound states (J'^ = 0~, 1~). In this case, the equation for 
the L- and N-component decouple, and one may verify that if eq.(^) possesses a 
pseudoscalar solution of the form L{q) = /(|q|), eq. (|7^) admits a vector solution cor- 
responding to the same orbital wave function multiplied by a constant polarization 
vector, N(p) = S'/(|p|), S = const.. (Here, we suppose that the 3-dimensional po- 
tential is a function of |p — q| only.) A similar degeneracy holds for arbitrary angular 
momentum and both parities, as is seen from the corresponding partial-wave equa- 
tions, eqs.( |A.9| , [A.IOD in appendix A.l. Thus, the bound state equations eqs.(^, |B^) 
naturally realize the heavy-quark spin symmetry in the limit rriQ oo. 

Let us briefly discuss the non-relativistic limit of this description of mesons. In 
this case we take in eq. (|6^) also i^jdpl) = and obtain =1,5^= 0. In addition, 
we approximate both quark energies by their non-relativistic values, £^„(|p|) = m„ + 



|p|V2m^ 



n 



,j. In this limit, eqs. (|7T|, [72D reduce to the well-known non-relativistic 



potential model for heavy quarkonia 

We have seen that the bilocal effective theory with instantaneous interaction pro- 
vides a unified description of both light and heavy fiavors. It should thus be well suited 



19 



to investigate the spectrum and decays of heavy-light bound states. Our intention here 
is not to perform an exhaustive calculation of the meson spectrum based on the Salpeter 
equations, eqs. (|63|, |6^) , but rather to obtain a quantitative estimate of the effects of 
dynamical chiral symmetry breaking on the heavy-light meson spectrum and decay 
constants. Specifically, we wish to demonstrate the phenomenological importance of 
describing chiral symmetry breaking self-consistently through the Schwinger-Dyson 
equation, eq. (|4^) , rather than by a constant mass for the light quark, eq.(^). To this 
end, we compare the masses and decay constants of the D and i?-mesons calculated 
using the momentum-dependent quark self energy from eq. f^i]) with those obtained 
using a constant light quark mass, eq. (^0|) . 

To describe the D and i?-mesons, we take as potential the sum of a Coulomb 
and a linear potential, eq.(|l6D, as has been used in the analysis of charmonium, cc 
1^. In our investigations here we use the same potential for D and S-mesons and 
ignore the running of a^. For the parameters we take the values a = 0.41 GeV and 
as = 0.39, which we determined by fitting the masses of the IS'- and 25'-state of 
cc {J/ip) as well as those of bb (T), using eqs.(|63|, |6^) in the non-relativistic limit, 
= 0,i?j(|p|) = _Ej(|p|) = nic + ||pp/mc. In this fit the quark masses were 



obtained as nic = 1.39 GeV and rrib = 4.79 GeV. 

The result of a calculation of the masses and decay constants of pseudoscalar D 
and 5-mesons and their first radially excited states with the above potential is shown 
in table ^ There, results are given obtained with both the dynamical self-energy for 
the light flavor and with a constant light quark mass. For the sake of comparison 
we have chosen the constant light quark mass equal to the dynamically generated 
quark mass, eq. (|49|) , at |p| = 0, which for this potential is m(0) = 0.082 GeV, cf. 
fig.3. We have also included values obtained with a typical light quark "constituent" 
mass, ruq = 0.33 GeV. As can be seen, the use of the constant quark mass instead 
of the dynamical self-energy leads to considerable changes in the meson masses and 
decay constants. The increase in the decay constant of the ground-state D and B 
meson can be explained by the sensitivity of F^^b to the higher-momentum part of 
the wave function, which is larger if a constant quark mass is used. This is seen also 
from fig. 7, which shows the radial wave functions in both cases. Note further that 
the D and i?-meson masses calculated with the dynamical quark mass agree rather 
well with the experimental values, much better than the ones obtained with either 
niq = const.. (If rriq is taken to be significantly smaller than the values in table 1, 
the bound state mass even increases.) However, one should not overemphasize the 
quantitative agreement, since we did not choose parameters such as to fit the strength 
of chiral symmetry breaking in the light sector. More observables (/,r, (qq)) should 
be involved in order to optimize the potential before drawing conclusions based on 
absolute values. Nevertheless, the results of table |l] show clearly that a dynamical 
description of chiral symmetry breaking is important, and that noticeable effects due 
to the dynamical nature of the quark masses occur already in the meson mass spectrum 
and decay constants. Our results confirm the conclusions of Kaburagi et ai, who 
investigated the dependence of heavy-light meson properties on the light quark mass 
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in the framework of the Dirac equation with a scalar confining potential |^ . It should 
be stressed, however, that the main point in having a description which takes in to 
account the dynamical breaking of chiral symmetry is a qualitative one. The clearest 
manifestations of the role of dynamical chiral symmetry in heavy-light systems may 
be seen not in the masses and decay constants but in more complicated observables, 
e.g. processes involving emission of pions. Such processes can be described using the 
effective bilocal meson action, eg . (^OD . 

Finally, we would like to comment on the ordering of the pseudoscalar meson decay 
constants. We find Fb to be slightly larger than Fd, for both the dynamical and 
constant light quark mass. This is in agreement with the calculation of Cea et al. 137 1, 



who use a relativistic potential model, but contrary to most non-relativistic results, 
cf. and references therein. In it is argued that relativisitc effects spoil the 

— 1/2 

simple proportionality Fh oc Mh ' . If we treat in our approach the heavy quark as 
non-relativistic, i.e., if we solve eqs.(|63|, [6^ ) with = and Eip = uiq + ||pp/mQ, 
we find Fd to be somewhat larger than Fb for both the dynamical self-energy and 
the constant quark mass for the light flavor. Thus, the pattern is reversed if the 
relativistic kinematics is abandoned. Note that our spectra were calculated using the 
same potential for D- and S-mesons; larger differences in the decay constants may 
occur with a running coupling constant. 



5 Summary and outlook 

In this paper we have presented a bilocal effective model, which describes in a unified 
way both light and heavy mesons. Chiral symmetry is broken spontaneously by the 
interactions which bind the mesons. The important feature of our model is that the 
potential kernel, eq.(|TDp, moves together with the bound state, which leads to a rela- 
tivistically covariant description of bound states. The bilocal meson action obtained 
from this model, eq. (pOD , provides the script for deriving equations for bound states and 
for the calculation of matrix elements. The correct relativistic kinematics is essential 
in describing meson decays and more complicated processes like semileptonic decays 



18, 23|. This covariant formulation, which constitutes a relativistic extension of the 



usual Coulomb gauge, should also be useful in the light quark sector, for example in 



describing the pion electromagnetic form factor |^ or processes like p — vrvr. 

The Salpeter equations for the meson wave functions can be simplified considerably 
by a Foldy-Wouthusen transformation, eq. (|^) . Of practical importance is the fact 
that this transformation does not lead to complications in the numerical solution of 
the bound state equations if the momentum-space Multhopp method is employed. The 
equations exhibit both chiral symmetry and heavy-quark spin symmetry in dependence 
on the chiral angle specifying the quark single-particle spectrum. 

With a phenomenological potential used in the description of charmonium, good 
results for the masses of the pseudoscalar D and i?-mesons are obtained if spontaneous 
chiral symmetry breaking is taken into account. Heavy-light meson masses and decay 
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constants are seen to be sensitive to the momentum dependence of the hght quark 
self-energy. This demonstrates the necessity of formuahing chiral symmetry breaking 
self-consistently and justifies the effort of solving the Schwinger-Dyson equation in the 
light flavor sector. 

Our aim here has been to set up the framework for describing heavy-light meson 
bound states in a bilocal effective theory. We have tried to demonstrate the viabil- 
ity of this scheme and to get a quantitative picture of the meson spectrum. Clearly, 
there is much room to improve and extend the results presented here. For example, 
an intermediate-range potential or transverse gluon exchange could be included in ad- 
dition to the Coulomb and confinement potential in order to increase the strength of 
chiral symmetry breaking in the light sector Furthermore, an interesting possi- 

bility would be to consider the effective meson action with the covariant instantaneous 
interaction in the heavy-quark limit and perform a long-wavelength expansion, as out- 
lined in p2|. This would allow one to make quantitative predictions for the lagrangian 
of heavy quark effective theory in dependence on the interaction potential assumed at 
quark level. Moreover, the covariant instantaneous formulation could be generalized 
to the particle-particle sector of the quark theory, eq.(^, to describe also baryons as 
bound states. 
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A The Salpeter equation 
A.l Partial wave decomposition 

In this section we decompose the Schrodinger-type equations in the rest frame, eqs.(|63|, 
|6^) , into equations for bound states of given total angular momentum and parity 
(c/. also [^). This is achieved by expanding the functions L{p), N(p) in spherical 



harmonics of total angular momentum, J. For the vector component it is convenient 
to use instead of the usual vector spherical harmonics the combinations 

YjM = pYjM = ofYjJ-lM — /5Yjj+iM, 

^JM = {J{J+l))-'/''Vf,YjM = (3Yjj_,M + aYjj+,M, (A.l) 

^JM = -t{J{J+l))-'/'pxVpYjM = YjjM, 

with 



a = J—^, (3 = \ a' + (3' = 1. (A.2) 

V2J + I' ^ V2j+r ^ ^ ' 

The new functions Y}^ have simple transformation properties under the operations 
p X Yjn^ and p-Y}^. They are orthogonal and normalized as the usual vector spherical 
harmonics. We thus write the component functions of eq. (|6^) as0 

L^'\p) = ^^F,m(p), N^^Hp) = E^^Y}^/(P) ik = l,2). (A.3) 
P t^i P 

For the angular matrix element of the potential kernel we use the definition 

^|rf^],|rffi,F;,,,,(p)v(p-q)FiM(q) = t;L(p,g)5ii'W- (a.4) 

Note that the matrix element is independent of M, and that VL{p,q) is a symmetric 
functions of the radial variables p, q. The general expressions for vl{p, q) corresponding 
to the power-like potentials used in eq.(O) are given in table El. From the formula 



PQ 
{2ny 



■ J dnpj dngY*,^,^,{-p)V{p-q)YjLM{cO = VL{p,q)5jj'5LL'SMM' (A.5) 

(L = J,J±1) 



we obtain the matrix elements 



(2vr) 



dQ,J dQ,Y)i,{p)vip - q)Y^M(q) 



a^vj-i + p^vj+i = vj{p,q) (A,p) = (l,l) 

+ = vj{p,q) (A,p) = (2,2) 

vj (A,p) = (3,3) 

a(3{vj.,-vj+i) = vj{p,q) (A,p) = (l,2),(2,l) 



(A.6) 



^In this section, p, q denote |p|, |q| 
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All other combinations vanish. 

With the above definitions it is straightforward to reduce eqs. (|63| , |6^ ) to the follow- 
ing two systems of radial equations, 

(2) (1) .00 . (1) (?) , 

M^J^ {p)-E/j^ {p) = dq[[c],c^Vj + sppjYj\q)+isppjn^j{q)} 

(?) Q) roc . ^ © ^ (?) . 

Mn^j (P) - EpTi^j (p) = dq{ [cjcpj + sjsjvj]nsj (g) - zsjs^ Vjij (q) | 

(A.7) 



and 



(1) (2) f (2) (2) 

Mn^j (p) - EpU^j (p) = dq[ [c^c^Vj + sp^Vj\n^j (g) + c^c^Vjn^j 

(1) (2) ( _ (2) I (2) 

Mn^j (p) - EpU^j (p) = dq[ [c^cjvj + spjvj]n2j (g) + cjc^ Vjn^j 



Here, vj = vj{p,q) etc.. The radial equations describe bound states of parity (— )'^^^ 
and (— )'^, respectively. 

In the heavy-quark limit, with the Foldy-Wouthuysen factors given by eg . (|69D , 
Cp = Cp, Sp = ±Sp, the systems of radial equations eqs.( |A.7| , [A.8| ) simplify to 

iM-E,(t)i.) ^ j:mc^c,[- -i:K,i)('' 

(A.9) 



n2j Jo \vjvj M vj r \n2j 



(A.IO) 



Here, we have dropped the indices on the wave function components in accordance 
with eg . (|70D . These reduced equations can be diagonalized by introducing new radial 
wave functions, 

aij = alj - iPn^j, bu = ariu + pn2j, 

a2j = - ioiHsj, 62J = -(3nij + an2j. 



For the transverse vector components this is eguivalent to using in eg.( A.ij| ) instead of 



the Y^jyj the ordinary vector spherical harmonics, Yjlm, as basis functions. Upon this 



the system egs.( |A.9| , |A.10|) becomes 



'---4-)'^' ^ r^'^'-K" -J'l::)''' 
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dq [cpCq I ^ ^^^^ \ + SpSq { 




(A.13) 



The new system eqs.( |A.l2 , A.13| ) exhibits the heavy-quark spin symmetry. For a given 



J the equations for ai^j+i and 02, j coincide with those for 62, j and respectively. 
This means that the bound states come in degenerate pairs with total angular momen- 
tum J and J + 1, for both positive and negative parity. For example, in this limit the 
0^ (pseudoscalar) meson is degenerate with the 1~ (vector), and the 0"*" (scalar) meson 
with the 1"*" (axial vector). 



A. 2 Numerical solution: the Multhopp method 



The radial equations, eqs. (|A.7| , |A.8|) , are in general singular integral equations with 



power-like or logarithmic singularities at g = j9. A simple and powerful numerical tech- 
nique to solve such equations directly in momentum space is the Multhopp method, 
which has been used frequently in the context of non-relativistic and relativistic con- 



stituent quark models |p4[. Of crucial importance for this method is the fact that for 
potentials with known singularities the integral of the basis functions with the singular 
part of the potential kernel can be split off and performed analytically, so that only 
finite integrals need to be evaluated numerically. In the radial equations eqs.( |A.7| , 
[A.8|) a new feature compared to the quark model is the presence of the momentum- 
dependent Foldy-Wouthuysen factors, c^,Sq , modifying the integration kernel. Here, 
we briefly show how the difficulties presented by these additional form factors can be 
circumvented and the singular part of the integrals be separated as usual. Thus, the 
Multhopp method can be used just as efficiently in the case of a non-trivial single- 
particle spectrum as in standard quark model calculations. 

For simplicity, we consider eq.( [A.7D for a = 0~ bound state; the generalization 
to the coupled equations for J > is straightforward. In this case Vq = 0, tJo = f 1 and 
eq.(|Xl) for (i^^^'^\p) = il)^'^\p) simplifies to 



(2) (1) .00 (1) 

Mi ' ip) - E,i ' ip) = dq [cJcJi;o(p, q) + sj3ji;i(p, g)]£ ' (g) (A.14) 

The Multhopp method consists in converting eq. ( |A. 14| ) into a matrix equation in mo- 
mentum space. The range of momenta (0, 00) is mapped onto the interval (0, tt) by the 
coordinate transformation 



p = Atan|6', q = Atan|x. (A. 15) 

where A is a scale parameter with dimensions of momentum. One then chooses as 
basis functions the finite set (2/7r)^/^ sinz^, i = 1, . . . N. The expansion of £^^'"^^9) is 
equivalent to interpolating i^^''^\9) at a set of angles 9k = kn / (A^ + 1), /c = 1, . . . A^, the 
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so-called Multhopp angles. By making use of the orthogonality relation for the basis 
functions one obtains from eg. ( [A. 141) a 2N x 2iV-eigenvalue equation for 



N (i) a\ (A 



^2- 

B 

k=l 



J2B.:t'ie,) = Mi' (9,), (A.16) 



with 



Q) 2 ^ (i) 

B^, = E{e,)6,, + j;^Y.^simd,B'{z,e,), (A.17) 

(1) \ 
B'it,e,) = / dx- ^[cfcpoie,x) + sfs^v,ie,x)]simx (A.18) 

Jo COS 2 "X^ 

Here, -^(6*), Cg", sj, v^iQ, x)-, x) are related to the corresponding functions 

in eq.( |A.14|) by the transformation, eq. ([A.15| ). The solution of eq.( |A.16D directly yields 



the momentum-space wave function at a discrete set of points [^. The calculations 
quoted in section 4 were performed with = 20 ... 30. 

The main input to eq. (|A.16|) are the Multhopp integrals of the potential kernel, 
eq.( |A.18D . In particular, the singularity of the integration kernel at g = p is absorbed 



in the integral over the basis functions, eq. ([A.18| ). We now show that by a simple 



rearrangement of terms the part of eq.( |A.18D containing the singularity can be sepa- 
rated even in the presence of the Foldy-Wouthuysen factors. The Coulomb potential 
produces only an integrable logarithmic singularity in the radial equation, c/. table |^, 
and poses no further problems. For the linear potential the singularity is of the type 

VL{p,q)im = 7 — ^— + log. div. terms. (A.19) 
[P - q) 

Note that the leading singularity is independent of the angular momentum, L. Let us 
consider the contribution of the leading singularity, eq.( [A.19D , to the radial integral of 



eq. (|A.14|) . As the leading singularities oivo{p, q) and q) are identical, the singular 



part of the integral in eq.( A.14| ) is given by the principal value integral 



{q-p) 



Since the radial wave functions satisfy £*-^'^^(0) = ^'^"'^'^-'(oo) = 0, eq. ([A.20|) can be 
integrated by parts, which gives 
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This can be rewritten as 



dq d Q, 



Jo {q — p) dq 



-i ' (q) (A.22) 



+ r {[44 + 44 - ii) + 144 - 44/^ (^) ^(-1. T -2,)} 

Of the integrals in eq.( |A.22| ) only the first one is singular, which does not involve 
the Foldy-Wouthuysen angle in the integrand, while the second one is non-singular. 
Thus, the radial integrals involving the linear potential can be rearranged such that the 
Foldy-Wouthuysen factors, which are in general known only numerically, enter only in 
non-singular integrals. In other words, the Foldy-Wouthuysen transformation does not 
modify the leading singularity of the radial equation, as is to be expected. After the 
substitution, eq. (|A.15|) , the Multhopp integral corresponding to the singular part of 
eq. (|A.22|) can be performed analytically as usual ||2^, because the Foldy-Wouthuysen 



factors do not enter. The remaining integrals of eq.( [A.22D , which involve the Foldy- 
Wouthuysen factors, lead to Multhopp integrals with at most logarithmic singularities, 
which can be calculated efficiently with standard Fourier transform routines. 

The above argument can easily be extended to the full equations eqs.( |A.7| , |A.8| ) for 



J > 0. In fact, it follows from eq.( |A.19| ) and eq.( [A.2D that the leading singularities of 



tJjip, q) and g), as defined in eq. ( |A.6|) , are equal to that of f g), and that the 
leading singularities cancel in vj{p,q). 



For a pure oscillator potential the entire Multhopp integral, eq.( [A.18D , can be eval- 



uated analytically. In this case the Multhopp method provides an alternative to solving 
the differential equation for the radial wave functions state by state [^, as the solution 
of eq.( |A.16| ) immediately gives the ground state and the excited state spectrum. 

The scale parameter A is chosen to be of the order of the r.m.s. 3-momentum of 
the bound state wave function. Eigenvalues and eigenfunctions should not be sensitive 
to this choice. We emphasize that this method is not a variational approach, i.e., no 
variation is performed with respect to A, in contrast to the procedure of pi, |35|. 



B Normalization and decay constants 

In this appendix we evaluate the matrix elements needed for the normalization of 
the bound state amplitude and for the calculation of the pseudoscalar meson decay 
constants within the bilocal field approach. We shall derive explicit expressions in 
terms of the bound state wave functions in the rest frame. 

For the normalization of the bound state amplitude we consider the matrix element 
of the free part of the quark loop, eq. (^7D , between on-shell bound states H ~ {qiqj), 

{H{V'^)\W^^\H{Vh)) = {27:)%V'^-VH){2uj'^2uJnVl'G-\VH\ (B.l) 
G-\Vh) = -i\NJ:i^{Vh)GiT{Vh)G2]. (B.2) 
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The trace implies integration over the loop momentum. We evaluate ^ ^ in the rest 
frame of the bound state, Vh = (Mj/,0). Writing T{q\VH) = r(q), r(q|P/^) = r(q| - 
Vh) = r(q), we have 



■tr 



r(q)G,(g - ip^)r(-q)G,(g + IVh] 



(2^ 



■tr 



_(2vr)4 
r(-q)A,+ (7or(q)7o)A, 
E,-Mh 



(B.3) 



+ 



r(-q)A,_(7or(q)7o)A, 



E„ + Mh 



The numerators here can be rewritten as 



tr[r(-q)A^±(7or(q)7o)Ay = 

-tr [(Sri(q)r(-q)Sri(q)) A± {S.\^)V{^)St\c^)) A 



(B.4) 



where we have used the relation between the free and the rotated projectors, eg. (^T]) . 
We now use the Salpeter equation, in the rest frame, eq. (|60D , and the corresponding 
equation for r(q) to express eq.( [B.3| ) in terms of the "undressed" wave function in the 
rest frame, eq. (|56|) . 



d^q 



E, - M^)tr [A- l'(-q) A+ ^(q)] (B.5) 



(27r)3 

+ {Eq + M^^)tr [A+ l'(-q) A- ^(q)]}- 



In a general frame, Mh is to be replaced by and the integration in eq. (p.5| ) goes 

over the transverse momentum, q^. The bound state amplitudes are normalized by 
the requirement that the dispersion relation be the same as that of a free relativistic 
particle. 



(B.6) 



Note that for a co-moving instantaneous interaction, eq.(0), the interaction part of 
the effective action, eq. (|20D , does not contribute to the normalization matrix element 
because depends on V only through the unit boost vector, r], and V^{d/dV^)rj = 0. 
From eq.( [B.5| ), the normalization condition in the rest frame can thus be expressed as 

Nr r d^v ( ~Q 00 oT oo 

1 = ^y^{tr[A-Vl/(-q)A+vI/(q)]-tr[A+vI;(-q)A-Vl/(q)].} (B.7) 



In terms of the components of the wave function, eq. 



this reads 



2Nr 



H 



d^p 
(2^ 



_^{l)*_^(2) _^ _^,(2)*_^(1) _^ n(1)*.n(2) _^ ]\[(2)*.]\j(2) 



(B.8) 
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where L^^)* = L(i)*(q) = ^(^^(-q), etc.. 

We now consider the pseudoscalar meson decay constants. The matrix element of 



eq.(|3^ for the decay of a meson H ~ (qiqj) into a leptonic pair is given by 



iHVLW^lmVH)) = {2i^f5^'\VH-VL)^{lu\l,\Q)T^{VH). (B.9) 

J^'iVu) = -iN,j^,ii,[O^GM-\VH)nq^\VH)GM+\VH)]. (B.IO) 

By integrating over the parallel component of the loop momentum, this matrix element 
can be evaluated in terms of the wave function of the "moving" bound state, eq. (|5^) . 



:F'[Vh) = N^Jj^,tT,[0^^iq^\VH)], (B.ll) 

where "^{q^lVn) = ^(^'"'"1 ~ 'Ph)- The decay constant is then read off by comparing 
eq.( |B.9| ) with the general definition 

MVlW^'LIH.A'Ph)) = (27r)V(P^-P,)^F.,n(^H/.|0). (B.12) 



Note that for 3-dimensional interactions of the form eq.(p!OD in general different values 
for the timelike and spacelike part of the pion decay constant are obtained [^. 

In this paper we only consider the decays of pseudoscalar mesons (JP = 0~) at rest, 
Vh = {Mh, 0). In this case the wave function has the form 



|(q) = {L«(q)-7o^('Hq)}f, (B.13) 

as can be seen from the partial-wave expansion in appendix A. Inserting this expression 
with eq. (|56|) into eq.( p.ll|) and calculating the trace we obtain the decay constant of 
a pseudoscalar meson H ~ (qiqj) at rest. 

Here, the wave functions in the rest frame L^^''^\q) satisfy eqs.(0, and are normal- 



ized according to eq. (|B.8|) . In particular, in the isospin limit, z/jg = = ^{^tt — (j){\q\)) 



the decay constant becomes 



For light mesons (vr, K) it is customary to define the decay constant as /tt = Fj,/ \/2, 
etc.. 
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Tables 





u. y iiQjiiiiLjCXi 


m — 082 


m — 33 




mo 


2.03 


2.33 


2.34 


1.87 


Fd 


0.16 


0.21 


0.25 


< 0.29 


rriD' 


2.54 


2.85 


2.89 




Fd' 


0.16 


0.16 


0.19 




rriB 


5.38 


5.64 


5.65 


5.28 


Fb 


0.18 


0.24 


0.29 




rriB' 


5.83 


6.11 


6.14 




Fb' 


0.21 


0.23 


0.26 





Table 1: The masses and decay constants of the 0~ D- and i?-mesons and their first 
radially excited states for a linear plus Coulomb potential with a = 0.41 GeV, = 
0.39, rric = 1.39 GeV and nib = 4.79 GeV. All energies in GeV. Masses and wave 
functions were determined from eg. (|03D , decay constants from eq.( B.14 ). Columns 
1-3 refer to different descriptions of the light quark spectrum. 1: dynamical chiral 
symmetry breaking, eq.(Q), with mP = 1 MeV, 2: constant light quark mass, eq. (|50D , 
mq " ' "'^ - ^ - 

experimental bound on Fd in column 4 is from ref. ^ . 



m(|q| = 0) = 0.082 GeV, 3: constant light quark mass, rriq = 0.33 GeV. The 
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V{r) 
1/r 



Vip - q) 

Air 



|p-q| 



Stt 



(27r)X^^'nP - q) 




VLip,q) 
-Ql{x) 

TV 

l]_dQL 
71 pq dx 

L{L + l 

9 



6{p-q) 



Table 2: The Fourier transforms and the angular matrix elements, eq. (|A.4|) , corre- 
sponding to the power-like potentials used in eg. (p!6|) . The linear potential involves a 
contact term, C = Svr / IqI""^, which cancels the IR-divergence. Here, Ql{x) is the 
Legendre function of the second kind, with x = [p"^ + q^)/2pq. 



Figure captions 

Fig.l: The chiral angle, 0(|p|), from the Schwinger-Dyson equation for a light flavor 
{wP = 1 MeV). Solid line: linear plus Coulomb potential (cr = 0.41 GeV, = 0.39), 
dotted line: pure oscillator potential (Vq = 0.247 GeV). The dashed line shows the 
value for a constant heavy quark mass, mq = rric = 1.39 GeV. 

Fig. 2: The quark energy, i^dpl), from the Schwinger-Dyson equation for a light flavor 
(m° = 1 MeV). The potentials are those of fig.l. The dashed line shows the energy for 
a constant heavy quark mass, mg = rric = 1.39 GeV. The straight solid line indicates 
the asymptotic behavior, E = |p|. 

Fig. 3: The dynamical quark mass, m(|p|) of eq. (|49D , for a light flavor (m° = 1 MeV). 
The potentials are those of fig.l. 

Fig. 4: The squared pion mass, m^, in the chiral limit, as a function of the current 



quark mass, 



m 



The potentials are to those of fig.l. Solid line: linear 



plus Coulomb potential, dotted line: oscillator potential. 

Fig. 5: The decay constants of the pion, /^r, and its first radially excited state, in 
the chiral limit. The potentials are those of fig.l. Solid lines: linear plus Coulomb 
potential, dotted lines: oscillator potential. Note the different scales for and /^r'. 
Here, U = F^/V2, f^i = F^,/V2. 



Fig. 6: The radial wave functions of the pion, i''^^'^-'(|p|), for a linear plus Coulomb 
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potential {a = 0.41 GeV, a, = 0.39) and m° = = 1 MeV. Solid lines: dotted 
lines: l^'^^ . Shown are the ground and first excited state. The single-particle spectrum 
is determined by the Schwinger-Dyson equation. 

Fig. 7: The (l)-component of the radial wave function of the D-meson, ^'^^^(Ipl), for a 
linear plus Coulomb potential (a = 0.41 GeV, as = 0.39) and rric = 1.39 GeV. Solid 
line: dynamical quark self-energy for the light flavor (m" = IMeV), dotted line: 
constant light quark mass, niq — 0.082 GeV. We have not shown it is very close 
to due to the large mass of the charmed quark. 
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